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A method and a device for detecting an 
abnormality of a heat exchanger, 
and the use of such a device 

5 

The present invention relates to a method and a 
device for detecting an abnormality of a heat ex- 
changer exchanging heat between a first fluid flow 
flowing in a conduit and a second fluid flow flowing 

10 along a flow path, said conduit and said flow path 
. each having an inlet and an outlet, and the use of 
such a device. 

Heat exchangers are an important part of many 
plants and systems, especially refrigeration or heat 

15 pump systems. These heat exchangers and their effi- 
ciency are of crucial importance in such systems, and 
it is therefore important to monitor functioning of 
the heat exchangers to be able to detect abnormality 
of the heat exchanger, so measures can be taken to 

2 0 remedy any defects. 

By abnormality of the heat exchanger is meant 
that the heat exchanger does not exchange as much en- 
ergy as expected, i.e. the fluids do not experience 
the cooling or heating they should. This may be due 

25 to fouling of heat exchanger, in that a layer of 
scale, dirt or grease is deposited on. the heat ex- 
changing surface or surfaces leading to reduced heat 
exchange, as this layer will usually act as an insu- 
lating layer. Another possibility is that there is 

30 insufficient fluid flow because of dirt or the like 
blocking or restricting flow through the heat ex- 
changer. Both situations lead to higher power con- 
sumption, because the system must work at a higher 
load than a system working with heat exchangers 

35 within the normal range. Further in the event of ad- 
verse working conditions with high heat exchange de- 
mand and a relatively small temperature difference 
between the fluids, it may be impossible to meet the 
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demand, which in some systems may have devastating 
effect . 

Often the abnormality will not be detected be- 
fore an adverse working condition is experienced, in 
5 that the demand cannot be met, e.g. leading to an in- 
crease of the temperature of a system, which should 
be kept at a specific temperature. An example of such 
a system is a refrigerated display cabinet in a shop, 
where strict legislation in most countries prescribes 

10 that when food is not kept below a maximum tempera- 
ture, it must be discarded, which of course is expen- 
sive and devastating for the business. Likewise large 
computer systems are often kept in air-conditioned 
rooms, as an excessive temperature may increase the 

15 risk of a computer crash, which may entail a high 
risk of data loss and lost man-hours. 

Common provisions for detecting abnormality of 
a be*t- exchanger inoludo baoic vioual inspection at 
regular intervals to check for dirt at the inlet of 

2 0 the heat exchanger. Often the heat exchangers are 

placed so inspection is difficult, and hence such in- 
spection is labour consuming. Further an abnormality 
may arise at different intervals and quite quickly, 
e.g. in the event of material blocking the inlet to 
25 the heat exchanger. This means that to provide a rea- 
sonable degree of security against heat exchanger ab- 
normality, it is necessary to inspect the heat ex- 
changers often. Further a visual inspection of the 
outside of the heat exchanger may not be effective in 

3 0 assessing whether the internal heat exchanging sur- 

faces are subject to fouling etc. causing a reduced 
heat exchange . 

Another known way of detecting abnormality of a 
heat exchanger is by direct flow measurement. A di- 
35 rect flow measurement requires delicate and expensive 
equipment, such as hot wire anemometers or the like, 
and a plurality of flow measurement devices should be 
used to gain useful information on the overall flow 
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field. It has also been proposed to assess the flow 
based on pressure sensors, but such pressure sensors 
are also expensive, and to gain useful information on 
the overall flow field a plurality of pressure sen- 
5 sors should be used. A further disadvantage of these 
methods are that they can only be used to establish 
if there is restricted flow in the heat exchanger, 
not the situation where the flow is normal, but heat 
exchange is reduced, e.g. because of fouling of the 
10 heat exchanger surface. 

It is an object of the .present invention to 
provide a method to enable early detection of an ab- 
normality of a heat exchanger. 

This object is met by a method comprising the 
15 steps of: 

establishing at least one parameter representa- 
tive of the temperature conditions of the heat ex- 
changer, 

establishing a second fluid inlet temperature, 
20- establishing <a parameter indicative of expected 

heat exchange between the heat exchanger and the sec- 
ond fluid, 

processing the heat exchanger temperature, the 
second fluid temperature and the parameter indicative 

25 of expected heat exchange for establishing an esti- 
mated second fluid outlet temperature, 

employing the estimated second fluid outlet 
temperature for evaluating the heat exchange between 
the first and second fluids by comparing the esti- 

3 0 mated second fluid outlet temperature, or a parameter 
derived therefrom, with a reference value. By this 
method is provided a convenient way of assessing the 
functioning of the heat exchanger primarily based on 
parameters of the first fluid, which means that a 

35 minimum of sensors are needed for providing informa- 
tion regarding the second fluid, and an automatic de- 
tection of heat exchanger abnormality is made possi- 
ble. Further this method makes it possible to detect 
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restricted flow as well as fouling of the heat ex- 
changer . 

According to an embodiment, the reference value 
is a predetermined second fluid outlet temperature. 
5 An even more reliable, alternative method can 

be obtained when the method comprises the steps of 
using the estimated second fluid outlet temperature 
for establishing a second heat rate of the second 
fluids for evaluating the energy balance of the second 

10 heat rate of the second fluid compared to a first 
heat rate of the first fluid, as 1 an evaluation based 
on the assumption of energy balance will take the in- 
fluence of further parameters into account. 

According to an embodiment, the method com- 

15 prises establishing the second rate of heat flow of 
the second fluid by establishing an estimate of a 
second fluid mass flow and a specific enthalpy change 
of the second fluid across the heat exchanger based 
on the estimated second fluid outlet temperature and 

2 0 the second fluid inlet temperature. 

According to an embodiment-, the method com- 
prises establishing the first rate of heat flow by 
establishing a first fluid mass flow and a specific 
enthalpy change of the first fluid across the heat 
25 exchanger based on parameters representative for 
first fluid inlet and outlet temperatures, and the 
condensation pressure. 

A direct evaluation of the heat exchange is 
possible, but may however be subject to some disad- 

3 0 vantages, e.g. because of fluctuations or variations 

of the parameters in the refrigeration or heat pump 
system, and according to an embodiment, the method 
comprises establishing a residual as difference be- 
tween, the first heat rate and the second heat rate. 
35 It may also be possible to evaluate the heat 

exchange by direct evaluation of the estimated outlet 
second fluid outlet temperature, but this may however 
be subject to some disadvantages, e.g. because of 
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fluctuations or variations of the parameters in the 
refrigeration or heat pump system, and according to 
an alternative embodiment, the method comprises es- 
tablishing a residual as difference between the esti- 
5 mated and predetermined second fluid outlet tempera- 
ture . 

To further reduce the sensibility to. fluctua- 
tions or variations of parameters in the system and 
be able to register a trend of heat exchanging, the 
10 method comprises providing an abnormality indicator 
by means of the residual, the abnormality indicator 
being provided according to the formula: 

_ + *P When S MJ-l + S M,i > 0 

where s Mi is calculated according to the follow- 
15 ing equation: 

where Ci is a proportionality constant, /i 0 a 
first sensibility value, and /i a second sensibility 
value . 

20 Another aspect of the invention regards a heat 

exchanger abnormality detection device for a heat ex- 
changer exchanging heat between a first fluid in a 
conduit and a second fluid in a flow path, where the 
device comprises a first estimator estimating a heat 

25 exchanger temperature, a first intermediate memory 
means storing the heat exchanger temperature, a tem- 
perature sensor measuring the second fluid inlet tem- 
perature, a second intermediate memory means storing 
the second fluid inlet temperature, a second estima- 

3 0 tor establishing a parameter indicative of expected 
heat exchange between the heat exchanger and the sec- 
ond fluid, a third intermediate memory means storing 
the parameter indicative of expected heat exchange, a 
processor establishing an estimated second fluid, out- 
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let temperature based on said heat exchanger tempera- 
ture, said second fluid inlet temperature, from the 
first and second intermediate memory means, respec- 
tively, and, from the third intermediate memory 
5 means, the parameter indicative of expected heat ex- 
change, and a comparator comparing the estimated sec- 
ond fluid outlet temperature, or a parameter estab- 
lished on basis thereof, with a reference value. 

An embodiment of the device further comprises 

10 memory means for storing at least one parameter from 
the processor, whereby a device is obtained which may 
operate on the basis of previously stored data. 

Although applicable to heat exchangers in gen- 
eral, it is found that the device is particularly 

15 suited for an embodiment, where the heat exchanger is 
part of a vapour-compression refrigeration or heat 
pump system comprising a compressor, a condenser, an 
expansion device, and an evaporator interconnected by 
conduits providing a flow circuit for the first 

20 fluid, said first fluid being a refrigerant. 

According to an embodiment, the heat exchanger 
is the condenser, which is particularly difficult to 
monitor, as the refrigerant in the condenser is pre- 
sent in three different phases, namely as superheated 

25 gas, a mixture of gas and liquid and sub cooled liq- 
uid. 

According to an embodiment, the second fluid is 
air, which is the most common type of second fluid 
. for refrigeration or heat pump systems as outlined 

30 above, and for which direct measurement of fluid pa- 
rameters involves some special problems. Further the 
air used normally is the ambient air, which may con- 
tain different kinds of pollution, which may deposit 
on the heat exchanger. 

35 Specifically the evaporator may be part of a 

refrigerated display cabinet positioned within a 
building and the condenser is positioned outside the 
building, which is a special example, where the de- 
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vice according to the invention may be of particular 
value . 

A third aspect regards use of a detection de- 
vice as outlined above, where the detection device is 
5 used for detecting fouling of the heat exchanger 
and/or detecting insufficient flow of the second 
fluid. 

In the following, the invention will be de- 
scribed in detail by way of example with reference to 
10 the drawing, where 

Fig. 1 is a sketch of a refrigeration system, 
Fig. 2 is a schematic sectional view of a heat 
exchanger , 

Fig. 3 is a schematic end view of the heat ex- 
15 changer, 

Fig. 4 is an example of a temperature profile 
in the heat exchanger, 

Fig. 5 is a schematic log p,h-diagram of a re- 
frigerant , 

2 0 ' Fig. 6 is a curve of estimated and measured 

outlet temperature of a condenser, 

Fig. 7 is a curve of a residual, 

Fig. 8 is a curve of an abnormality indicator, 

and 

25 Fig. 9 is an enlarged portion of the curve ac- 

cording to Fig. 8. 

In the following reference will be made to a 
heat exchanger in a simple refrigeration system, al- 
though the principle is equally applicable to a heat 

3 0 exchanger in other heat exchanging systems, and as 

understood by the skilled person, the invention is in 
no way restricted to a refrigeration system. 

In Fig. 1 is shown a simple refrigeration sys- 
tem 1 comprising a compressor 2, a condenser 3, an 
35 expansion device 4 and an evaporator 5, which are 
connected by a conduit 6 in which a refrigerant is 
circulating. In vapour -compress ion refrigeration or 
heat pump systems the refrigerant circulates in the 
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system and undergoes phase change and pressure 
change. In the system 1 a refrigerant gas is com- 
pressed in the compressor 2 to achieve a high pres- 
sure refrigerant gas, the refrigerant gas is fed to 
5 the condenser 3 (heat exchanger) , where the refriger- 
ant gas is cooled and condensates, so the refrigerant 
is in liquid state at the exit from the condenser 3, 
expanding the refrigerant in the expansion device 4 
to a low pressure and evaporating the refrigerant in 

10 the evaporator 5 (heat exchanger) to achieve a low 
pressure refrigerant gas, which can be fed to the 
compressor 2 to continue the process. 

A specific example for the heat exchanger is a 
condenser 3 of a refrigeration system for a frozen 

15 food storage cabinet, or a refrigerated display cabi- 
net for shops. In Fig. 2 can be seen a schematic sec- 
tional view of a cross-flow heat exchanger with a 
first fluid flow 7 in a conduit 6 and a second fluid 
flow 8 in a flow path 9, Given the example of a con- 

2 0 denser 3 of a refrigeration system 1 as mentioned 
above, the first fluid is a refrigerant and the sec- 
ond fluid would normally be air. The refrigerant en- 
ters the condenser 3 as a superheated gas, which dur- 
ing the passage of the condenser 3 is cooled by the 

25 air flowing around and past the hot conduits 6 con- 
taining the refrigerant, so the refrigerant gas is 
cooled to condensation temperature, condensates and 
leaves the condenser 3 as sub cooled liquid. To ob- 
tain and maintain a flow of air through the condenser 

30 3, the condenser 3 is normally provided with a fan 
(not shown) , which can be running constantly, in an 
on-off mode or with a varying speed. 

Typically a condenser 3 of such a system is 
placed outside the shop, often on the roof, because 

35 if it is placed inside, it would lead to temperature 
increase in the shop, and normally the outside tem- 
perature is lower than the inside temperature. How- 
ever placing the condenser 3 outside has the drawback 
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the condenser 3 may be exposed to clogging up or 
fouling because of dirt, grease, leaves, newspapers, 
etc restricting the air flow 8 or reducing heat 
transfer from the refrigerant to the air, and further 
5 the heat exchanger is difficult to reach and to in- 
spect. The intervals of clogging up may be very ir- 
regular due to weather conditions, e.g. precipita- 
tion, wind direction etc., pollution, seasonal 
changes, such as leaf fall, which makes it difficult 

10 to provide proper inspection intervals. 

Fig. 3 • is a schematic end view of a heat ex- 
changer, which in the given example is the condenser 
3, the inlet of which is partly covered by leaves. 
This means that the air flow through the flow path 9 

15 of the condenser 3 is restricted, and hence that heat 
exchange is reduced. To be able to detect such an ab- 
normality, an air outlet temperature estimator is 
proposed using an air flow dependent thermo conduc- 
tivity a. The value of the thermo conductivity is de- 

2 0 .pendent on the given heat exchanger and can be estab- 

lished at start up of the heat exchanger. It is found 
that the value of this parameter is not critical, and 
the value may be established based on empirical val- 
ues or values supplied by the manufacturer of the 
25 heat exchanger. The thermo conductivity is flow de- 
pendent, and for a heat exchanger comprising a fan 
forcing air through the heat exchanger, the thermo 
conductivity a can be expressed as 

a = a 0 co°' s , (1) 

3 0 where <x> is the speed of the fan, and oc 0 is the 

thermo conductivity at no-flow condition. 

An estimate for the air temperature through the 
condenser 3 can be defined using the thermo conduc- 
tivity a, which for a constant airflow gives 

35 dt air (y) = a(W^O0-'*G0)*' (2.) 

where y denotes the distance from the air 
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inlet, and y is represented as a normalized parame- 
ter, i.e. is the distance relative to the total 
length of the flow path, so the outlet is at 1, 
tcond t surf(y) ^ S the surface temperature of the condenser 
5 heat exchanging surface at position y, and t air (y) is 
the air temperature at position y. t coadsUff could be es- 
tablished by direct measurement using a temperature 
sensor. Such sensors, however, are expensive and es- 
pecially when placed outside they are subject to er- 

10 rors . It is hence preferred to establish an estimate 
of the surface temperature of the condenser based on 
evaluation of parameters of the refrigerant. 

In the condenser the refrigerant is present in 
three different phases: in a region at the refriger- 

15 ant inlet, the refrigerant is in gas phase and more 
or less superheated, in another region, in which the 
refrigerant condensates at a constant temperature, 
the refrigerant is present as a mixture of gas and 
liquid, and in a third region, the refrigerant is 

20 liquid and more or less sub-cooled. 

In the following, t condsur/ (y) is in the two-phase 

and liquid regions assumed to be equal to the conden- 
sation temperature of the refrigerant. For the gas 
phase region ^^ tJur/ (v) is assumed to be the mean of the 

25 refrigerant gas temperature and the condensation tem- 
perature . 

For the region of the heat exchanger, where the 
refrigerant is present as sub-cooled liquid, which is 
the region from y=0 to y=yi, it is assumed that the 

30 temperature increases linearly with a gradient k ir 
which is found to produce an adequately accurate tem- 
perature profile for most purposes. Estimating as 
mentioned that the condenser surface temperature 
equals the condensation temperature, the temperature 

35. increment can be found as: 

dt air O) = a{t liquid + k x y - t air (y))dy ( 3 ) 
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where k x is a constant describing the tempera- 
ture gradient in the sub-cooled region, and t liquid is the 
refrigerant temperature at the refrigerant outlet. 

For the second, two-phase region from y=yi to 
5 y=Y2/ estimating as mentioned that the condenser sur- 
face temperature equals the condensation temperature 

( t cond , surf = t cond ) / 

dt+(y) = a(t^-t+(y)}dy (4) 
For the superheated gas region, i.e. the region 
10 from y=y 2 to y=y3=l, the temperature is assumed to 
vary linearly with a gradient k 2 , and the equation is 
" estimated as 

dt aif (y) = a(t cond +k 2 (y-y 2 )-t air (y))dy (5) 
where y 2 denotes the end of the two -phase re- 
15 gion, and k 2 denotes a mean temperature gradient for 
the superheated gas phase. 

To obtain the air temperature in the condenser, 
the above equations (3) , (4) and (5) are integrated, 
and hence: 

2 0 For the sub -cooled region (0<y<yi) 

where t a ir(0) is the air temperature at the 
inlet to the flow path 9, i.e. the ambient tempera- 
ture . 

25 For the two-phase region (yi<y<y 2 ) 

tM=t^ y* c * _, ° (7) 

, For the superheated gas region (y2<y<l (=y3 ) ) 

It is hence possible to estimate an air outlet 

3 0 temperature using these equations. Parameters needed 

for the air outlet temperature estimate are the air 
inlet temperature, the temperature of the refrigerant 
at the inlet and outlet, the condensation temperature 
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of the refrigerant, estimates for y x and y 2 , ki and k 2 . 
It is found that for many condensers approximately 5% 
of the heat exchange is in the first region where the 
refrigerant is present as sub-cooled liquid, approxi- 
5 mately 75% of the heat exchange takes place in the 
second region, i.e. the part of the condenser where 
the refrigerant is changing phase from gas to liquid, 
and the remaining approximately 20% of the heat ex- 
change takes place in the region of the condenser, 

10 where the refrigerant is present as superheated gas. 
The value of k x can be established more or less em- 
pirically based on y lf t liquid and t cmd , whereas k 2 can be 
established more or less empirically based on y 2/ the 
refrigerant outlet temperature, and the overall 

15 length of the flow path. The air outlet temperature 
can thus be obtained primarily based on parameters of 
the refrigerant, and these refrigerant parameters 
will normally already be known, as most modern re- 
frigeration systems comprise a controller of the re- 

20 f j_iycj-aliuii system with sensors constantly measuring 
these parameters. Fig. 4 illustrates an example of 
the temperature profile of the cross -flow condenser. 
If a less accurate answer is satisfactory, it is pos- 
sible to use a simplified mod'fel taking e.g. only the 

25 two-phase region and the superheated gas region into 
account, or even only the two-phase region, where 
most of the heat exchange takes place. 

The estimated air outlet temperature can then 
be compared with a measured air outlet temperature 

3 0 obtained by a temperature sensor at the air outlet. 
When the heat exchanger experiences an abnormality, a 
significant estimation error occurs, which can be 
used to trigger an alarm signal. 

Although this approach of comparing the air 

3 5 outlet temperature . directly with a measured tempera- 
ture may be convenient and adequate in some systems, 
a more stable and reliable result . can be obtained 
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when basing the evaluation on the assumption of en- 
ergy balance of the heat exchanger. However, a direct 
outlet temperature is rarely convenient and moreover 
temperature sensor measuring air outlet temperature 
5 will seldom be present, so there is a need for an al- 
ternative approach. 

The energy balance of the condenser can be 
stated as: 

10 where Q Air is the heat taken up by the air per 

time unit, i.e. the rate of heat flow delivered to 
the air, and Q Ref the heat removed from the refriger- 
ant per time unit, i.e. the rate of heat flow deliv- 
ered by the refrigerant. 

15 . The basis for establishing the rate of heat 

flow of the refrigerant (Q Ref ) i.e. the heat delivered 

by the refrigerant per time unit is the following 
equation: 

Q,Ref = ™Ref ( k Ref, in ~ Kef, out ) ( 1 0 J 

20 where m Ref is the refrigerant mass flow. h Refout 

is the specific enthalpy of the refrigerant at the 
condenser outlet, and h Refin is the specific enthalpy 
of the refrigerant at the condenser inlet. The spe- 
cific enthalpy of a refrigerant is a material and 

25 state property of the refrigerant, and the specific 
enthalpy can be determined. The refrigerant manufac- 
turer provides a log p, h-diagram of the type accord- 
ing to Fig. 5 for the refrigerant, wherein the ther- 
modynamic cycle of a refrigeration system is sketched 

3 0 for illustration. From I to II, the refrigerant gas 
is compressed in a compressor, from II to III, the 
refrigerant is cooled in a condenser from a state of 
superheated gas " to condensation and further to a 
state of sub-cooled liquid. From III to IV, the re- 

3 5 frigerant is expanded in an expansion device to a 
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lower pressure, where the refrigerant is present as a 
mixture of liquid and gas. From IV to I, the refrig- 
erant is heated in an evaporator so at point I at the 
entry to the compressor, the refrigerant is com- 
5 pletely gaseous. 

With the aid of this diagram the specific en- 
thalpy difference across the condenser can be estab- 
lished. For example to establish h Ref in with the aid of 

a log p, h-diagram, it is only necessary to know the 

10 temperature and the pressure of the refrigerant at 
the condenser inlet (T Re f,m and Pcond/ respectively) . 
Those parameters may be measured with the aid of a 
temperature sensor and a pressure sensor. 

Similarly, to establish the specific enthalpy 

15 at the condenser outlet, two measurement values are 
needed: the refrigerant temperature at condenser out- 
let (T Re f,out) and the pressure at the condenser outlet 
(Pcond) / which can be measured with a temperature sen- 
sor and a pressure sensor, respectively. 

20 Instead of the log p, h-diagram, it is natu- 

rally also possible to use values from a chart or ta- 
ble, which simplifies calculation with the aid of a 
processor. Frequently the refrigerant manufacturers 
also provide equations of state for the refrigerant, 

25 so a direct calculation can be made. 

The mass flow of the refrigerant may be estab- 
lished by assuming solely liquid phase refrigerant at 
the expansion device entry. In refrigeration systems 
having an electronically controlled expansion valve, 

3 0 e.g. using pulse width modulation, it is possible to. 
determine the theoretical refrigerant mass flow based 
on the opening passage and/or the opening period of 
the valve, when the difference of absolute pressure 
across the valve and the subcooling (T v ,i n ) at the ex- 

35 pansion valve entry is known. Similarly the refriger- 
ant mass flow, can be established in refrigeration 
systems using an expansion device having a well-known 



WO 2004/036170 



PCT/DK2003/00070J 



15 

opening passage e.g. fixed orifice or a capillary- 
tube. In most systems the above-mentioned parameters 
are already known, as pressure sensors are present, 
which measure the pressure in the condenser 3 . In 
5 many cases the subcooling is approximately constant, 
small and possible to estimate, and therefore does 
not need to be measured. The refrigerant mass flow 
through the expansion valve can then be calculated by 
means of a valve characteristic, the pressure differ- 

10 eritial, the subcooling and the valve opening passage 
and/or valve opening period. With many pulse width 
modulated expansion valves it is found for constant 
subcooling that the theoretical refrigerant mass flow 
is approximately proportional to the difference be- 

15 tween the absolute pressures before and after and the 
opening period of the valve. In this case the theo- 
retical mass flow can be calculated according to the 
following equation: 

<ef=^<P cond -P £vap yOP (ID 

20 where is the absolute pressure xn the don- 

denser, P Evap the pressure in the evaporator, OP the 
opening period and a proportionality constant,, 

which depend on the valve and subcooling. In some 
cases the subcooling of the refrigerant is so large, 

2 5 that it is necessary to measure the subcooling, as 
the refrigerant flow through the expansion valve is 
influenced by the subcooling. In a lot of cases it is 
however only necessary to establish the absolute 
pressure before and after the valve and the opening 

30 passage and/or opening period of the valve, as the 
subcooling is a small and fairly constant value, and 
subcooling can then be taken into consideration in a 
valve characteristic or a proportionality constant. 
The value of the mass flow is not critical, and an- 

35 other possibility is to establish the mass flow from 
the compressor directly based on empirical values 
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e.g. data supplied by the manufacturer of the com- 
pressor and the absolute pressure before and after 
the compressor. 

Similarly the rate of heat flow heat of the air 

5 (Qxir)' i>e - the heat taken up by the air per time 
unit may be established according to the equation: 

QAir=™Air(Kir,ou t - h Air,J ^ 

where m Mr is the mass flow of air per time 
unit, h Airin is the specific enthalpy of the air before 
10 the condenser, and h Airout is the specific enthalpy of 

the air after the condenser. 

The specific enthalpy of the air can be calcu- 
lated based on the following equation: 

h Ai =l.QQ6-t+x(250l + l.&-t\[h] = kJ/kg (13) 
15 where t is the temperature of the air, i.e. 

T a ir,in before the condenser and T a i r ,out after the con- 
denser, x denotes the absolute humidity of the air. 
The absolute humidity of the air can be calculated by 
the following equation: 

20 x = 0.62198 ^ (14) 

PAmb-Pw 

Here p w is the partial pressure of the water' 
vapour in the air, and p Amb is the air pressure. p Amb 
can either be measured or a standard atmosphere pres- 
sure can simply be used. The deviation of the real 

25 pressure from the standard atmosphere pressure is not 
of significant importance in the calculation of the 
amount of heat per time unit delivered by the air. 
The partial pressure of the water vapour is deter- 
mined by means of the relative humidity of the air 

3 0 and the saturated water vapour pressure and can be 
calculated by means of the following equation: 

p w =p WJSat -XH (15) 
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Here RH is the relative humidity of the air and 
p WSat the saturated pressure of the water vapour. p WSat 

is solely dependent on the temperature, and can be 
found in thermodynamic reference books. The relative 
5 humidity of the air can be measured or a typical 
value can be used in the calculation. 

When equations (10) and (12) is set to be 
equal, as implied in equation (9) , the following is 
found : 

10 Air (Kir, Out " Kir, In ) = ™Ref ( h R ef, In " K ef, Out ) ( 1 6 > 

From this the air mass flow. m Air can be found by 
isolating m Mr : 

(^Ref, In ^^Ref, Out) / 1 -7 \ 

171 Air = m Ref * 1 7 ' 

(J 1 Air, Out ~k Air Jn ) 

Assuming faultless air flow this equation can 
15 be used to evaluate the operation of the system. In 
many cases it is recommended to register the air mass 
flow in the system. As an example this air mass flow 
can be registered as an average over a certain time 
period, in which the refrigeration system is running 
20 under stabile and faultless operating conditions. 
Such a time period could as an example be 100 min- 
utes. This estimated air mass flow found as an aver- 
age under stabile and faultless operating conditions 

is denoted m Air . 

25 A certain difficulty lies in the fact that the 

signals from the different sensors (thermometers > 
pressure sensors) are subject to significant varia- 
tion. These variations can be in opposite phase, so a 
signal for the estimated air outlet temperature or 

30 the energy balance is achieved, which provides cer- 
tain difficulties in the analysis. These variations 
or fluctuations are a result of the dynamic condi- 
tions in the refrigeration system. It is therefore 
advantageous regularly, e.g. once per minute, to es- 

35 tablish a value, which in the following will be de- 
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noted "residual", based on the energy balance accord- 
ing to equation (9) : 

r = Q Air -Q«* (18) 
so based on the equations (10) and (12) , the 
5 residual can be found as: 

r =W2 Air (^Air, Out ~~ ^Air, In )~™Ref (^R flf , In ~ ^R ef, Out ) ( 1 9 ) 

where m Air is the estimated air mass flow, which 
is established as mentioned above, i.e. as an average 
during a period of faultless operation. Another pos- 

10 sibility is to assume that rh Air is a constant value, 
which could be established in the very simple example 
of a condenser having a constantly running fan. Even 
in systems with ^variable flow capacity, such as sys- 
tems having a plurality of fans, which can be acti- 

15 vated independently, or in systems incorporating one 
or more fans running with variable speed, e.g. using 
a frequency converter, a fair estimate of the mass 
flow can be established. The estimated mass flow can 
be found by establishing the number of currently con- 
nedted rans, i.e. how many fans are connected, and/yx 
the speed of the fans, to thereby establish the flow 
capacity of the connected fans, e.g. by use of em- 
pirical values. 

The estimated air outlet temperature can simi- 

25 larly by evaluated by providing the residual as the 
difference between the estimated air outlet tempera- 
ture and a predetermined air outlet temperature. The 
predetermined air outlet temperature may be measured 
directly or may be obtained as an empirical value. 

30 In a refrigeration system operating fault- 

lessly, the residual r has an average value of zero, 
although it is subject to considerable variations. To 
be able to early detect a fault, which shows as a 
trend in the residual, it is presumed that the regis - 

3 5 tered value for the residual r is subject to a Gaus- 
sian distribution about an average value and inde- 
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pendent whether the refrigeration system is working 
faultless or a fault has arisen. 

In principle the residual should be zero no 
matter whether a fault is present in the system or 
5 not, as the principle of conservation of energy or 
energy balance of course is eternal. When it is not 
the case in the above equations, it is because the 
prerequisite for, the use of the equations used is not 
fulfilled in the event of a fault in the system. 

10 In the event of fouling of the condenser sur- 

face, the thermo conductivity changes, so that a be- 
comes several times smaller. This is not taken into 
account in the calculation, so the estimated rate of 
heat flow of the air Q Air used in the equations "is 

15 significantly bigger than in reality. For the rate of 
heat flow of the refrigerant (O^), the calculation 
is correct (or assumed correct) , which means that the 
calculated value for the rate of heat flow of the re- 
frigerant ( Qkzf ) across the heat exchanger equals the 

20 rate of heat flow of the refrigerant in reality. The 
consequence is that the average of the residual be- 
comes positive in the event of fouling of the con- 
denser surface . 

In the event of a fault causing reduced air 

25 flow through the condenser (a defect fan or e.g. dirt 
covering the air inlet of the heat exchanger) the 
mass flow of air is less than the estimated value of 

the mass flow of air m Air used in the calculations. 

This means that the rate of heat flow of the air used 

30 in the calculations is larger than the actual rate of 
heat flow of the air in reality, i.e. less heat per 
unit time is removed by the air than expected. The 
consequence (assuming correct rate of heat flow of 
the refrigerant) , is that the residual becomes posi- 

35 tive in the event of a fault causing reduced air flow 
across the condenser. 
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To filter the residual signal for any fluctua- 
tions and oscillations statistical operations are 
performed by investigating the residual. 

The investigation is performed by calculating 
5 an abnormality indicator according to the following 
equation: 

{0,M>henS^+s Mti <0 (20) 

where s Mi is calculated according to the follow- 
. ing equation: 

^=4 r '- i T £ ) 1211 

where d is a proportionality constant, fi Q a 
first sensibility value, fi a second sensibility 
value, which is positive. 

In equation (2 0) it is naturally presupposed 

15 that the abnormality indicator S Mjt i.e. at the first 
point in time, is set to zero. For a later point in 
time is used according to equation (21) , and the 

sum of this value and the abnormality indicator. S § at 
a previous point in time is computed. When this sum 

20 is larger than zero, the abnormality indicator is set 
to this new value. When this sum equals or is less 
than zero, the abnormality indicator is set to zero. 
In the simplest case /x 0 is set to zero. \i is a chosen 
value, ' which e.g. establish that a fault has arisen. 

25 The parameter \i is a criterion for how often it is 
accepted to have a false alarm regarding heat ex- 
changer abnormality detection. 

When for example a fault occurs in that the air 
inlet of the condenser is covered by e.g. leaves, 

30 then the abnormality indicator will grow, as the pe- 
riodically registered values of the j , in average is 
larger than zero. When the abnormality indicator 
reaches a predetermined value an alarm is activated, 



WO 2004/036170 



PCT/DK2003/000701 



21 

the alarm showing that the air mass flow is reduced. 
If a larger value of /x is chosen, fewer false alarms 
are experienced, but there exist a risk of reducing 
sensitivity for detection of a fault. 
5 The principle of operation of the filtering ac- 

cording to equation (20) and (21) shall be illus- 
trated by means of Figs. 7 and 8, where the filtering 
is used on the residual found using energy balance, 
i.e. based on equation (18). In Fig. 7 the time in 

10 minutes is on the x-axis and on the y-axis the resid- 
ual r. Fig. 7 illustrates the emerging of a fault in 
that the condenser of a shop was subject to a sudden 
fouling at approximately t=2900 minutes. However, as 
can be seen the signal is subject to quite signifi- 

15 cant fluctuations and variations, which makes evalua- 
tion difficult, and the presence of a problem is 
really not evident before approximately t=5500 min- 

In Fig. 8, which represent the filtering of the 

2 0 data of Fig. 7 with means of the abnormality indica- 

tor according to equation (20), the time in minutes 
in on the x-axis and on the y-axis the abnormality 
indicator S. As can be seen, the heat exchanger was 
working properly until approximately t=2900 minutes, 
25 when a sudden fouling took place, and the abnormality 
indicator S rose. This is easier to see in Fig. 9, 
which is an enlarged portion of Fig. 8. In Fig. 9 the 
abnormality at approximately t=2900 minutes can be 
easily detected using the abnormality indicator S 

3 0 compared to using the residual or the air outlet tem- 

perature . 

A further advantage of the device is that it 
may be retrofitted to any refrigeration or heat pump 
system without any major intervention in the refrig- 
35 eration system. The device uses signals from sensors, 
which are normally already present in the refrigera- 
tion system, or sensors, which can be retrofitted at 
a very low price. 



